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Abstract 



We present the results of a high statistics lattice calculation of hadronic form factors 
C^ \ relevant for D— and i?— meson semi-leptonic decays into light pseudoscalar and vector 

mesons. The results have been obtained by averaging over 170 gauge field configurations, 
generated in the quenched approximation, at /3 = 6.0, on a 18^ x 64 lattice, using the 
0(a)-improved SW-Clover action. From the study of the matrix element < K~\J^\D^ >, 
we obtain /+(0) = 0.78 it 0.08 and from the matrix element < K*^\J fj\D^ > we obtain 
F(0) = 1.08±0.22, ^i(O) = 0.67±0.11 and ^2(0) = 0.49±0.34. We also obtain the ratios 
V{<S)/Ai{id) = 1.6 ± 0.3 and ^2(0)/^i(0) = 0.7 ± 0.4. Our predictions for the different 
form factors are in good agreement with the experimental data, although, in the case of 
^2(0), the errors are still too large to draw any firm conclusion. 

With the help of the Heavy Quark Eff'ective Theory (HQET) we have also extrapolated the 
lattice results to S-meson decays. The form factors follow a behaviour compatible with 
the HQET predictions. Our results are in agreement with a previous lattice calculation, 
performed at /3 = 6.4, using the standard Wilson action. 
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1 Introduction 

There is increasing evidence that quantitative calculations of weak decay amplitudes can 
be obtained by lattice QCD simulations. Over the last few years, semi-leptonic decays 
of heavy mesons have been studied on the lattice [01-@1- Among these, Z)— meson decays 
provide a good test of the lattice method, since the relevant CKM matrix element is well 
constrained by unitarity in the Standard Model, Vcs — 0.975. The main advantage of the 
lattice technique is that it is based on first principles only and it does not contain free 
parameters besides the quark masses and the value of the lattice spacing, both of which 
are fixed by hadron spectroscopy. Moreover, statistical and systematic errors in lattice 
simulations can be systematically reduced with increasing computer resources. 

In this work, we present a high-statistics study of pseudoscalar-pseudoscalar and pseu- 
doscalar-vector semi-leptonic form factors, performed on the 6.4 Gflops APE machine, by 
using the 0(a)-improved Clover Action ||8|, |9[, at /? = 6.0, corresponding to an inverse 
lattice spacing a~^ ~ 2 GeV. Our prediction for the form factors, which govern the 
D ^ K and D -^ K* amplitudes are given in the abstract and in table |3|. The central 
values are in remarkable agreement with the experimental results. However, in spite of the 
large statistics, our predictions still suffer, in some cases, from large statistical errors. A 
possible explanation for the errors' size is the use of a "thinning" procedure, that will be 
discussed in detail in sec. |3[ This procedure was adopted because our computer memory 
is not sufficient to store the necessary quark propagators. 

Following the suggestion of ref. , we have also tried to extrapolate the form factors to 
B ^ IT, p decays, using the scaling laws predicted by the Heavy Quark Effective Theory 
(HQET) [^]. The final results have large uncertainties, because the statistical errors 
amplify in the extrapolation. It is reassuring that the present results are in good agreement 
with those obtained, with the same method, in ref. ||5|, by using the standard fermion 
Wilson action, at a smaller value of the lattice spacing, corresponding to /? = 6.4. This 
gives us confidence on the feasibility of the extrapolation. It should be noticed however, 
that the form factors are obtained near q"^ = q'^^x > where q^ is the momentum of the lepton 
pair. In order to predict the form factors at all q'^, on current lattices and in the range 
of the heavy quark masses explored so far, we can only assume pole dominance (or any 
other simple g^-dependence). This strongly biases the final results, e.g. the determination 
of the form factors at q^ = 0. In this respect, the situation is not very different from the 
quark model approach. In order to improve the accuracy of the predictions, it is necessary 
to be able to work with heavier quark masses and to increase the range of q'^. This can 
only be achieved by going to larger lattices. The approach followed here indicates that 
5— meson semi-leptonic decays will be a fruitful area of lattice investigations in the near 
future. 



2 Computation Details 

Semi-leptonic decays {D -^ K,K*; B -^ D,D*; B ^ tt,p) are described in terms of six 
independent dimensionless form factors, four of which are important for decay rates into 
light leptons (see e.g. Q). For each form factor, the relevant information can be expressed 
in terms of its value at q^ = and its q^ dependence. In this section several details of our 
calculation are given. 



Figure 1: The effective mass, defined as M{t)a = ln(C(t)/C(t + a)) as a function oft/a. Mpy 
and Mpy refer to the heavy-light or the light-light pseudoscalar and vector mesons respectively. 
C{t) is the generic, zero momentum two-point correlation function. 

We have obtained our results, by averaging over 170 gauge field configurations, gener- 
ated at /3 = 6.0, on a volume 18^ x 64. On this set of configurations, we have computed the 
quark propagators by using the 0(a)-improved SW-Clover action |^, ^. We have consid- 
ered 3 different values of the Wilson hopping parameter for the light quarks, Ki = 0.1425, 
0.1432 and 0.1440, and 4 different values for the heavy quarks, Kr = 0.1150, 0.1200, 
0.1250 and 0.1330. The values of K, corresponding to the zero quark mass and to the 
strange quark mass (obtained by fixing the pion and the kaon pseudoscalar masses) are 
Kcr = 0.14545(1) and Kg = 0.1435(2) respectively. The value of the charm quark mass, 
obtained from the Z)— meson mass, corresponds to Kch = 0.1219(17). The inverse lattice 
spacing, obtained by using the mass of the />— meson to set the scale, is a~^ = (1.95 it 0.08) 
GeV. All two-point functions have been fitted in the interval t/a = 12 — 28 and 15 — 28 
for light-light and heavy-light mesons respectively. In fig. |l], we give examples of the 
effective mass, as a function of time, in the two cases. The statistical errors have been 
estimated with the jacknife method, by decimating 10 configurations from the total set 
of 170. Preliminary results obtained with a smaller set of gauge field configurations have 



been reported in ref. |11|. 



To extract the form factors, we have computed two- and three-point correlation func- 
tions and followed the procedure explained in ref. P|. The matrix elements have been 
computed using the pseudoscalar density as source for D and K mesons and the local 
vector current for the K* mesonQ. For the weak current we have used the lattice im- 



^Throughout this paper, K , K* and D are conventional names to denote the light pseudoscalar and vector 



proved local vector and axial vector currents, according to the definition of refs. |12, ^]. 
The corresponding renormalization constants, Zy and Z^, have been fixed to the values 
Zy = 0.824 and Z^ = 1.06, as determined non-perturbatively, using light quark correla- 
tion functions [O, llH]. This choice will be justified below. 

All matrix elements have been computed by inserting the L'-meson source at a time 
distance {t£) — tK,K*)/(i = 28. The position of the light meson is fixed at the origin and 
we have varied the time position of the weak current in the time interval tj/a = 10 — 14. 

To study the (^^-dependence of the form factors, we have computed the matrix elements 
in two different kinematical configurations: in the first case we have taken the D— meson 
at rest, i.e. po = 0, &nA pk,K'' = 27^ /La- (0,0,0), (1,0,0), (1,1,0), (1,1,1) and (2,0,0); 
in the second case we have chosen pr, = 2t:/ La ■ (1,0,0) and pk,k* = 2iT/La • (0,0,0), 
(1,0,0), (—1,0,0), (0,1,0). We have therefore nine independent momenta. We have also 
computed other correlation functions, which are equivalent under the cubic symmetry. 
All such equivalent cases have been averaged together. 

In order to extract the current matrix elements from the three-point functions, we have 
used two different procedures, denoted by "analytic" and "ratio" methods, discussed in 
detail in ref. P]. Within the statistical errors, the above methods are expected to agree, 
up to 0(a) effects. In the present calculation, we find that the two methods yield slightly 
different results, the differences varying between 2% and 8%. We have taken into account 
these differences in the evaluation of the final error (see below). 

In order to obtain the form factors at q"^ = 0, for quark masses corresponding to the 
physical D (B) and K (vr) or K* (p) mesons, we have extrapolated the form factors, both 
in masses and momenta. The following procedure has been used: 

1) At fixed heavy quark mass and meson momenta, Pk,k*, the generic form factor F 
(F = f^, Ai, A2 and V) has been extrapolated linearly in the light quark mass, to 
values corresponding to the strange [D -^ K,K*) or massless (D -^ vr,p) quarks. 

2) F has been extrapolated in the mass of the heavy meson, to the D— and S— meson 
masses, using the dependence expected in the HQET (see eqs.(|^ below). In order 
to evaluate the stability of these results with respect to a different extrapolation, 
we have also extrapolated F in the mass of the heavy meson Mp according to the 
"naive" expression F = a + (3/Mp. The differences between the two methods are 
discussed later on. 

3) In order to obtain the form factors at g^ = 0, we have only used the points with 
Pd = and pk,k* = 2iT/La ■ (1,0,0), which have been extrapolated by assuming 
meson pole dominance, F{q^) = F(0)/(1 — q^ /Mf). This reduces the uncertainty 
of the extrapolation because, in most of cases, the point with p^ = 27r/La • (0, 0, 0) 
and pk,k* = 27r/La • (1,0,0) corresponds to the smallest q^. Mt is the mass of 
the lightest meson exchanged in the t-channel. Thus, the vector and axial meson 
masses have been used for the extrapolation of f~^,V and ^1,^2 respectively. Mt 
has been computed on the lattice, over the same gauge field configurations, at the 
same heavy and light quark masses used for the three-point functions. In order to 
obtain the physical D— and 5— meson masses, we have fitted the vector-pseudoscalar 
mass difference AM = Mp* — Mp (and similarly for the axial and scalar cases) as 
AM = Am + Bpi/Mp. The results, extrapolated to the strange and massless light 



mesons and the heavy meson respectively. 



Md* 


Md** 1++ 


Md** 0++ 


Md,,* 


Mb,** 1++ 


Mb,** 0++ 


1.95 ±0.01 


2.29 ±0.34 


2.03 ±0.19 


2.05 ±0.01 


2.40 ±0.34 


2.13 ±0.19 


Mb. 


Mb** 1++ 


Mb** 0++ 


Mb,* 


Mb," 1++ 


Mb,"0++ 


5.32 ±0.01 


5.99 ±0.62 


5.46 ± 0.32 


5.46 ±0.01 


6.13 ±0.62 


5.60 ±0.32 



Table 1: Masses (in GeV) of the vector, axial and scalar excitations for the D and B mesons 
as determined from our simulation. The experimental pseudoscalar masses Mb and Mb are 
used as an input. These masses have been used to extrapolate the form factors to q^ = 0. 

quarks, are reported in table ||. These values have been used to extrapolate the 
form factors at g^ = 0. In this extrapolation, the precise value of Mj is relatively 
unimportant. For example we have verified that, by using the vector meson mass in 
all cases, the results change by only a few per cent. 

The g^-dependence of l/f~^{q'^) and l/Ai^q"^) is compared with the meson dominance 
predictions in fig. ^. The values of the inverse form factors are given as a function of the 
dimensionless variable q'^/Mf, for the values of the heavy and light hopping parameter, 
Kh = 0.1250 and Ki = 0.1432. The lines in the figures represent the pole dominance 
expectations, with the pole masses computed on the lattice. In order to reduce the error 
due to the extrapolation in g^, the values of the form factors at zero momentum transfer 
have been obtained, by fitting only the point closest to q^ = 0, as explained above. 
The other points give, however, important information on the g^-dependence of the form 
factors. As shown in the figure, the q'^ behaviour of the two form factors is compatible 
with pole dominance predictions. Similar conclusions can be reached also for the form 
factors l^(g^) and ^2(9^)- Notice that the axial form factors determined by using QCD 
sum rules |l^, |l^ do to follow the pole dominance behaviour. 

3 Main Sources of Errors 

In this section, we briefly describe the main sources of error which are present in our 
calculation, besides the quenched approximation. 



3.1 The "Thinning" Procedure 

"Thinning" means that, when computing the correlation functions, one uses only one 
point out of Nth, in each spatial direction. In our case Nth = 3. Thus, for example, the 
Fourier transform of the two-point correlation function is defined as 



C{t,p)=NrY.^-'^'''^C{t,XT) 



(1) 



XT 



where xt = 0,3, 6, ... in each direction and A/r is a suitable normalization factor. This 
procedure is necessary when, as in our case, the computer memory is not sufficient to store 
the full quark propagators. There is a systematic error introduced by thinning, because we 
cannot eliminate high momentum components in the correlation functions. For Nth = 3, 
it is possible to show that, for each spatial direction, two higher momentum components 



Figure 2: l//_|_(g^) and 1/Ai(g^) as a function of the dimensionless variable q^/M^. The heavy 
and light Wilson parameters correspond to Kh = 0.1250 and Ki = 0.1432 respectively. The 
lines represent the pole dominance approximation. 

(pi and P2), besides the smallest one (po), give a contribution to the correlation function, 
pi and p2 are related to po by the simple relation: 

Pm = Po + -^ • ( -^ I -"i , m = l,2 (2) 

The systematic error introduced by thinning is expected to be negligible at large time 
distances, since the contribution of the unwanted higher energy states is exponentially 
suppressed in time. On the other hand, the resulting signal may be noisier, because we 
use a small sample of the lattice points. 

In the case of the two-point correlation functions, we have been able to directly compare 
"thinned" and "non-thinned" correlation functions, computed on the same set of gauge 
configurations, and for momenta p = 2ti/ La ■ (0,0,0) and (1,0,0). For the pseudoscalar 
and vector meson correlations, no observable statistical or systematic effect has been 
detected. However, the statistical noise introduced by thinning may be larger, when we 
compute three-point correlation functions, simply because, in the latter case, we have 
to thin twice. Thus, the thinning procedure could be responsible for the quite large 
statistical errors which have been found, in spite of the high statistical sample used in 
this calculation. Another reason could be the small spatial volume of our lattice. 



3.2 0{g^a) and 0{a^) Effects 

Using the Clover action, discretization errors are of 0{g^a) and 0{a?'). In the case of hght 
quarks, these effects have been shown to be much smaher than in the case of the standard 
Wilson action ||T^. In the charm quark mass region, at /? ~ 6.0, however, discretization 
errors may still be important also in the Clover case. 

An estimate of the lattice artefacts can be obtained by comparing values of lattice 
renormalization constants, computed from different matrix elements, in a non-perturbative 
way [^] . In ref . |l^ , the renormalization constants of the local vector and axial vector 
currents, V^ = V'i7^V'2 and A^ = '0i7^75'02; have been determined, at /3 = 6.0, as a 
function of the quark masses, mi and m2, using the Ward identity method p^ , p!^ . A 
variation of about 10 — 15%, on the values of Zy and Z^, has been observed, for heavy 
quark masses between m^a = 0.3 and 0.9, which is the range used in the present study. 



The renormalization of the quark fields, proposed in ref. |1£] in order to reduce 0{mHo) 
effects, enters, in our case, only at Oirn^jO?), since we are using an improved action and 
improved operators & 0]. Since the effects that we have observed are linear in niHa, 
the KLM correcting factor |jl^ is of no help in our case. We interpret the residual mass 
dependence as an effect of 0{asa). Similar results have been observed by the UKQCD 
Collaboration, which works with a Clover action, at /? = 6.2 Q. For the above reasons, we 
have used, through this paper, the non-perturbative values of Zy and Za, as determined 
in refs. |14, ^], for light quark masses. A 10 — 15% of systematic error is then certainly 



present in our final results, due to residual 0{asa) effects. 

3.3 Extrapolations of the Form Factors 

As explained in the previous section, the form factors, computed in a range of heavy 
quark masses around the charm, must be extrapolated in 1/Mp and g^, in order to obtain 
predictions for the physical form factors at zero momentum transfer. There is also an 
extrapolation in the light quark mass, but this is quite smooth and unlikely to be a source 
of an important uncertainty, within the present statistical accuracy. The extrapolation of 
the form factors in 1/Mp is more delicate. There are arguments, based on HQET, which 
allow the expansion of the form factors in inverse powers of the heavy meson mass Mp. 
On the basis of HQET, up to 0(1/Mp) and up to logarithmic corrections, one expects 
the following behaviour for the relevant form factors [|1C| ]: 

i = ^ , X 1-1 — = ^,^ X 1 -I 



(3) 
^i^/Mp = 7i X fl + -^) ^^=72xfl + — 



Mp) y/M^ ' V Mp, 

The expansions given in eqs. (^) are valid, in the limit of large heavy quark mass, at 
fixed momentum p of the light meson (in the frame where the heavy meson is at rest) 
when \p\ <C Mp. The above conditions are always satisfied for q^ax^ when the initial and 
final mesons are both at rest. In our simulation, they are also satisfied for the points 
corresponding to p = 211 /La -(1,0, 0) (these have been used in order to obtain all our final 
predictions). 

In sec. Q, we show that the dependence of the form factors on the heavy meson mass, 
Mp, is compatible with the HQET predictions. In order to evaluate the stability of the 



results with respect to a different extrapolation, we have also used a "naive" scaling law 
of the form: 



P+ .. . . _„/. , /? 



r,V,Ai,A2 = ax^l + j^j (4) 

For D— mesons, we find that the differences between the different extrapolations, eq.(^) or 
(^, are quite small (< 2%) and completely negligible with respect to the statistical errors. 
The results for S— meson decays will be reported in sec. ^. In this case the differences 
are larger, of the order of 10 — 20%, but still small with respect to the statistical errors, 
which amplify in the extrapolation. Thus, we are not able to distinguish between the two 
different behaviours. 

We now turn the discussion to the q'^ extrapolation. In the range of heavy quark 
masses considered in this study, the g^-dependence of the form factors is compatible with 
pole dominance predictions. This range corresponds approximately to a heavy meson 
mass of the order of the D— meson mass and we also have points at q^ ~ 0. Therefore, 
as far as the D decays are concerned, the extrapolation to g^ = does not represent 
an important source of theoretical uncertainty. The extrapolation in 1/Mp to the B- 
meson mass results in momentum transfers close to the maximum one. In order to then 
predict the form factors at small q^, we have been forced to assume pole dominance. 
Since the range of the extrapolation is very large^ this strongly biases our final results, 
e.g. the determination of the form factors at q'^ = 0. In particular, we have used pole 
dominance also in those cases, where it is in contradiction with the results of QCD sum 



rules calculations [16, 17]. In order to improve the situation, it is necessary to be able to 
work with heavier quark masses and to increase the range of q^. This can only be achieved 
by going to larger lattices (at larger values of /?). 

3.4 The "Analytic" and "Ratio" Methods 

In order to extract the weak current matrix elements from the three-point functions, we 
used two different procedures, denoted by "analytic" and "ratio" methods. The two pro- 
cedures are discussed in detail in ref. (H. The "ratio" methods means that, at each fixed 
time distance, we divide the three-point correlation function by the two relevant two-point 
functions of the D and K{K*) meson (with corresponding momentum) in order to cancel 
the exponential time-dependence of the three-point correlation function. The "analytic" 
method differs from the previous one because, instead of dividing by the two-point correla- 
tion functions, computed at different momenta, we divide by the corresponding analytical 
expressions, with the source matrix elements and the meson energies (computed from the 
meson masses, see eqs.(Q) and (^) below) taken from the fit of the two-point functions, at 
zero momentum. 

When both the initial and final mesons are at rest these two methods are practically 
equivalent and lead to almost identical results. However, when the meson momenta are 
different from zero, within the statistical fluctuations, the two methods are expected to 
agree only up to 0[a) effects. 

At 13 = 6.2 and with a lattice volume 24^ x 48, it was found that discretization 
errors seem to be reduced, by using, for the two-point correlation functions, the "lattice" 



^ While the typical Vy o'^ o^i'" lattice is at most 1.6 GeV, the extrapolation in the heavy quark mass will 
bring us to y/q^ = 4.5 GeV {q^ = Ml - 2Mba/M2 + (2 7r/18a)2 + M^), and similarly for B ^ p. 





/+(0) 


^(0) 


Ai(0) 


analytic 


0.81 ±0.07 


1.07 ±0.21 


0.66 ±0.11 


ratio 


0.75 ±0.06 


1.09 ±0.22 


0.67 ±0.11 




^2(0) 


y(o)Mi(o) 


A2(0)Mi(0) 


analytic 


0.44 ±0.34 


1.63 ±0.29 


0.67 ±0.43 


ratio 


0.52 ±0.29 


1.62 ±0.29 


0.79 ±0.35 



Table 2: Semi-leptonic form factors at zero momentum transfer for D ^ K and D —>■ K* 
decays. The results have been obtained by using the method called "analytic" or "ratio" to 
extract the form factors. 



dispersion relation of a free boson 



C{t,p) 



z 



-Et 



2sinh^ 



where 



E = — arcsinh 
a 



sinh\^) + Y. sinY-^) 



=1,3 



2 



(5) 



(6) 



and p is the momentum of the meson [^0|. The same is true in our case. We have 
verified this point, by studying the ratio R{t) = C{t,p)/C{t,p). C{t,p) is the two-point 
correlation function of a meson with momentum p, as computed in our simulation. C{t,p) 
is the expression in eq. @, with Z and m taken from the fit of the zero momentum 
correlation to C{t,p = 0). At large time distances, we find that R{t) is close to one 
(typically 1.05 ± 0.01). This compares favourably to the value 1.10 ± 0.02 in the case 
where, instead of eq. (^), we use the standard expression 



C{t,p-) = ^e-^\ 



with 



E = Jm^ + \p\' 



(7) 



(8) 



For this reason, in all our analysis, we have fitted the two point functions to C{t,p), as 
defined in eq. (|5|). 

Using C{t,p), the "ratio" and "analytic" methods yield slightly different results. In 
order to exhibit the differences, we give in table || the corresponding sets of values for the 
form factors /+(0), V{0), ^i(O) and ^2(0), as well as for the ratios V/Ai and A2/A1 at 
zero momentum transfer, for the D ^ K and D -^ K* decays. From the entries in the 
table one can see that the differences are larger for the decay in the pseudoscalar channel 
and they are always between 2% and ~ 8%pi In our final results, the differences between 
the two methods have been added in quadrature with the statistical errors. Had we used 
the standard expression C{t,p), we would had found larger differences between the "ratio" 
and "analytic" methods, which in the worst case have been evaluated to be 11%, see also 
ref. m. 



■^ For the form factor A2, the difference is of the order of 15%. In this case, however, the statistical errors 
are so large that the possibility of drawing any conclusion is prevented. 



Reference 


f+(o) 


V(0) 


Ai(0) 


EXP 


Average [2T| 


0.77 ±0.04 


1.16 ±0.16 


0.61 ±0.05 


LAT 


This work 


0.78 ± 0.08 


1.08 ±0.22 


0.67 ±0.11 




LMMS 


PI 


0.63 ± 0.08 


0.86 ±0.10 


0.53 ±0.03 




BKS |4| 


0.90 ±0.08 ±0.21 


1.43 ±0.45 ±0.49 


0.83 ±0.14 ±0.28 




ELC H 

UKQCD § 

LANL 


0.60 ±0.15 ±0.07 
n fi7+o-07 

"J-O'-O.OS 

0.73 ± 0.05 


0.86 ±0.24 
1 m+O-30 

-'-■UJ-_o.l3 

1.27 ±0.08 


0.64 ±0.16 
0.66 ±0.03 


SR 


BBD |1| 


0.60+^:}^ 


1.10 ±0.25 


0.50 ±0.15 


QM 


WSB |2| 


0.76 


1.23 


0.88 




ISGW |23] 


0.8 


1.1 


0.8 




GS ll 


0.69 


1.5 


0.73 


Reference 


A2(0) 


V(0)/Ai(0) 


A2(0)/Ai(0) 


EXP 


Average [pl| 


0.45 ± 0.09 


1.90 ±0.25 


0.74 ±0.15 


LAT 


This work 


0.49 ± 0.34 


1.6 ±0.3 


0.7 ±0.4 




LMMS 


31 


0.19 ±0.21 


1.6 ±0.2 


0.4 ±0.4 




BKS P 




0.59 ±0.14 ±0.24 


1.99 ±0.22 ±0.33 


0.7 ±0.16 ±0.17 




ELC 1 
UKQCD 
LANL 


i 
PI 


0.40 ±0.28 ±0.04 

0.66t°:}° 

0.44 ±0.16 


1.3 ±0.2 

1 4+0-5 

1.83 ±0.09 


0.6 ±0.3 

0.9 ±0.2 

0.74 ±0.19 


SR 


BBD |T6| 


0.60 ±0.15 


2.2 ±0.2 


1.2 ±0.2 


QM 


WSB |22| 


1.15 


1.4 


1.3 




ISGW |23] 


0.8 


1.4 


1.0 




GS [m 


0.55 


2.0 


0.8 



Table 3: Semi-leptonic form factors for D ^ K and D -^ K* decays. The label EXP refers 
to the experimental results and the labels LAT, QM and SR correspond to lattice, quark model 
and sum rules calculations respectively. 

4 Physics Results 



4.1 D— Meson Decays 

Our best estimates for the form factors and partial widths are those given in the abstract 
and in tables y and §. In these tables, we present our results, together with other cal- 
culations and experimental determinations of the form factors. Our central values for 
the four relevant form factors are in good agreement with the experimental data (typ- 
ically to within less than 10%), although we still suffer, in some cases, from sizeable 
statistical errors, especially in the case of A2. Considering that lattice calculations do 
not contain free parameters, we find this agreement remarkable. We also observe (tables 
P and 0) that predictions from QCD sum rules calculations are in agreement with lat- 
tice calculations and experimental determinations, whereas quark models fail to describe 
the D — > K* decay. Besides the results reported in the tables, we have also obtained 
/+(0)//f (0) = 1.02 ±0.03. 



EXP 



LAT 



EXP 



LAT 



SR 



QM 



Reference 



Average |^, ^ 



This work 

LMMS i 

ELCi 
UKQCD i 




Reference 



Average M 



This work 

LMMS |] 

ELCi 
UKQCD i 



Ball llj 



WSB H 
ISGW M 



T{D ^ K) 



9.0 ±0.5 



9.1 ±2.0 

5.8 ±L5 

5.4 ±3.0 ±1.4 

7.0 ±1.6 ±0.4 



6.4 ±1.4 



8.5 
7.1 



T{D 



w 



vr 



0.9 



-0.2 



0.8 ±0.2 

0.5 ±0.2 

0.5 ±0.3 ±0.1 

0.52 ±0.18 ±0.04 



0.39 ±0.08 



0.72 
0.38 



T{D -^ K* 



5.1 ±0.5 



6.9 ±1.8 

5.0 ±0.9 

6.4 ±2.8 

6.0^ 



-0.8 
-1.6 



3.8 ±1.5 



9.7 
9.2 



T{D 



TT 



0.6 ±0.2 

0.40 ± 0.09 

0.6 ±0.3 ±0.1 

0.43 ±0.11 



0.12 ±0.03 



0.68 
0.46 



Tl{D ^ K*} 

Tt{D -^ K*) 

1.15 ±0.17 



1.2 ±0.3 

1.51 ±0.27 

1.4 ±0.3 

1.06 ±0.16 ±0.02 



0.86 ±0.06 



0.92 
1.09 



r(D. 



4.0 ±0.6 



6.4 ±1.1 
4.4 ±0.6 



7.9 



Table 4: Semi-leptonic partial widths (in units oflO^^s ^) for D ^ K, K* , it, p and (p, using 
Vcs = 0.975 and Vcd = 0.222. The ratio of the longitudinal to transverse polarization partial 



widths for D — >■ K* is also given. The experimental values ofT{D 



TT 



and T(Ds 



have been computed by taking the corresponding branching ratios and meson life-times from ref 
[26]. 
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Figure 3: Extrapolation of the form factors to the D- and B-meson masses (diamonds) using 
the predictions of the HQET. The points corresponding to the lightest heavy quark mass have 
not been used in the fits. 

4.2 Extrapolation to 5— Meson Decay 

At the values of lattice spacing currently used in numerical simulations, we are unable 
to study directly the b quarks. However, as discussed above, in order to obtain indirect 
information on i?— meson semi-leptonic decays, we can follow the strategy suggested in 
ref. H]: we study the form factors in the region of the charm quark mass and then 
extrapolate the results to the bottom mass by using the scaling behaviour predicted by 
the HQET, eqs. (y). In order to reduce the uncertainty due to the extrapolation, one 
could also compute the form factors in the static limit, i.e. the limit in which the heavy 
quark mass is infinite. This determination is not available yet. 

Our results show that the dependence of the form factors on the heavy meson mass, 
Mp, is compatible with the HQET predictions, see eq.(^), as first observed in ref. |^]. 
This dependence is shown in fig. |^, where the four relevant form factors, extrapolated 
to the chiral limit in the light quark mass, are given, as a function of 1/Mp (crosses). 
The form factors were computed with p£) = and pk,k* = 27r/La(l,0,0). The values 
interpolated/extrapolated to the D and B meson masses (diamonds) are also given. 

From the values of fig. y, we can compute the form factors at q'^ = 0. They have 
been obtained by assuming the pole meson dominance behaviour, using the meson masses 
given in table y. The corresponding results are presented in table y, labelled "b" , together 
with the results of other theoretical determinations. As discussed earlier, in order to 
evaluate the stability of these results with respect to a different extrapolation, we also 
give the values (labelled as "a") obtained with the naive scaling laws, eq. (^. In the 
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LAT 



Reference 

This work "a" 
This work "b" 

ELC "a" PI 
ELC "b" ll 



SR 



QM 



Ball 131 



WSB HI 
ISGW gl 



Reference 



LAT 



SR 



QM 



This work "a" 
This work "b" 

ELC "a" 1 1 
ELC "b" 1 1 



Ball 



WSB HI 
ISGW gg] 



f+(0) 

0.29 ±0.06 

0.35 ±0.08 

0.26 ±0.12 ±0.04 

0.30 ±0.14 ±0.05 



0.26 ±0.02 



0.33 
0.09 



A2(0) 



0.24 ±0.56 

0.27 ±0.80 

0.38 ±0.18 ±0.04 

0.49 ±0.21 ±0.05 



0.4 ±0.2 



0.28 
0.02 



V(0) 
0.45 ± 0.22 
0.53 ±0.31 
0.34 ±0.10 
0.37 ±0.11 



0.6 ±0.2 



0.33 

0.27 



V(0)/Ai(0) 



2.0 ±0.9 
2.6 ±1.9 
1.4 ±0.2 
1.6 ±0.3 



1.2 

5.4 



Ai(0) 
0.29 ±0.16 
0.24 ±0.12 
0.25 ±0.06 
0.22 ±0.05 



0.5 ±0.1 



0.28 
0.05 



A2(0)/Ai(0) 



0.8 ±1.5 

1.0 ±3.1 

1.5 ±0.5 ±0.2 

2.3 ±0.7 ±0.2 



1.0 
0.4 



Table 5: Semi-leptonic form factors for B ^ tt and p. The label "a" refers to the naive 
extrapolation in 1/Mp, eq. (3), and label "b" to the extrapolation given in eq. (4)- To extrapolate 
to zero momentum transfer we have used the masses of table 1. We have taken the ISGW form 
factors, as extrapolated to q^ = 0, in ref. [17]. 

case of the £)— meson the difference between the HQET scaling laws and naive scaling 
is immaterial. For the i3— meson, the differences are smaller than the statistical errors 
(typically 10 — 20%), so that we are not able to distinguish between the two behaviours. 
We find reassuring, however, that the results for i?— decays are in good agreement with 
the results of ref. Q, where the standard fermion Wilson action and a smaller lattice 
spacing, corresponding to /5 = 6.4, were used. Finally, we observe that we find values for 
Ai and A2 smaller than those obtained by QCD sum rules, cf. table ^. The reason is 
probably due to the fact that we have assumed pole dominance, whereas QCD sum rules 
find that the axial form factors are flat in q^. 

From the values of the form factor /+(0) given in table 0, and by assuming meson 
dominance for the g^-dependence, we can give an estimate of the B ^ it decay rate. To 
evaluate the errors, we have allowed the form factor to vary in all possible ways by one 
a within the statistical errors and to vary in all possible ways among the values obtained 
with different extrapolations in 1/Mp, i.e. flts "a" and "b". We obtain 



r{B -^ -Klui) = \Vub\ (8 ± 4) X lO^^s 



(9) 



that corresponds to the branching ratio 

B{B ^^lui) = \Vub\'' {12 ±Q) (10) 

when the value tb = (1.49 ± 0.12) x 10^^^ sec. is used for the 5— meson lifetime. On the 
other hand, with our present accuracy, the errors on the form factors for the i? — > /? decay 
are still too large to get an estimate of the corresponding branching ratio. A accurate 
determination of this quantity can only be achieved by working with heavier quark masses 
and by going to larger lattices. 
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